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Refinements to the bow and stern seal subroutines of the 
XR-3 Loads and Motion computer program are made. New lift 
fan performance maps are developed and presented.' Comparison 
studies of simulated performance of the original XR-3 model 
with the old and new fan maps are made under varying 
sea conditions. Measured roll and pitch transient 
behavior of the XR-3 test craft in calm water is compared 
with the simulated response of both the original and the new 
model. In addition, simulated performances of the original 
and new models in irregular ahead seas are compared. 
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I. INTRODUCTION 



A. BACKGROUND 

The speed limitation of the conventional displacement 
vessel has been recognized for many years. In an effort to 
effect a quantum increase in surface vessel speed the United 
States Navy has been developing various craft whose principal 
means of support is other than hydrostatic lift. 

One type of craft currently receiving much attention is 
the surface effect ship or craft. In this category are the 
air cushion vehicles and captured air bubble craft (CAB). 

The general nature of these craft and their construction is 
well presented by Robert L. Trlllo in Ref. 1. In these 
craft either all or a major portion of its support is derived 
from a pressure differential between the atmosphere and a 
plenum chamber v/hich is open at the bottom. The speed 
advantages gained by this type of vessel are derived from 
two principal characteristics. It does not waste energy 
displacing a large volume of v;ater and the frictional forces 
betvjeen the hull and water are greatly reduced by the craft 
riding on a cushion of air v/ith only a minimum of the craft 
structure actually in contact \vlth the water. 

Vessels of this type are generally separated into two 
categories — "Air Cushion Vehicles" v;hose weight is entirely 
supported by the pressure differential in the plenum chamber 
and "Captured Air Bubble Craft" whose v;eight is partially 
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supported by sidewalls which extend Into the water. It is 
the second type of vessel in which the United States Navy 
has become primarily interested for ocean going applications. 

The term Captured Air Bubble (CAB) is possibly a misnomer 
since in this type of craft the air does leak out and has to 
be continually replenished by supply fans. The leakage 
though is relatively slow when it is compared to the air 
cushion vehicle which has continuous air gap around the 
periphery of the craft and hence the fans that replenish the 
air in the plenumi chamber have to be larger and more powerful 
than those of a CAB of similar size. 

This thesis is concerned with one aspect of the develop- 
ment of the Captured Air Bubble vehicle and that is the 
sim.ulatlon of such craft utilizing a digital computer, 
specifically the Loads and Motion Program developed by 
Oceanics, Inc. 

The basic rigid body analysis and spatial relationships 
of the Loads and Motion program are well documented in 
Ref. 2. Also covered are the principal dynamic and static 
approximations used in developing the equations of motion 
for the craft in its six degrees of freedom. 

A copy of this program was delivered to the Naval 
Postgraduate School in October of 1972 which had been 
specifically tailored for the Bell (lOOB) 100 ton craft. 

Leo and Boncal in Ref. 3 succeeded in converting this program 
to represent the XR-3 test craft. 
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B. OBJECTIVE 



The purpose of this thesis was to take a second more 
detailed look into all of the aspects of seal modeling and 
the fan performance map with the idea of possibly making 
Improvements in these representations before testing the 
computer model of the XR-3 under dynamic conditions. 
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II. FAN PERFORMANCE MAP ANALYSIS AND COMPARISON 



A . INTRODUCTION 

As the surface effect ship is predominately dependent 
on the plenum lifting force for its performance, the 
characteristics of the lifting fans are important if valid 
studies of the craft in question are to be made. 

Leo and Boncal reported in Ref. 3 that the steady state 
operating conditions of the XR-3 were obtained by increasing 
the leakage area of the craft until the proper equilibrium 
pressure was maintained in the main plenum. This area was 
twenty-six square feet. As the cross-sectional area of the 
plenum is less than twenty square feet this leakage area was 
out of proportion to properly represent the leakage area of 
the XR-3. 

While operating above hump speed the XR-3 test craft has 
been observed to have a gap underneath the sternseal of 
from four to five inches. Slight leakage does occur at 
the bowseal as reported by Ref. 3 but the total leakage area 
is estimated to be not over four square feet. 

To appreciate the effect this leakage area has on computed 
craft performance an understanding is required of the method 
by v/hich the pressure is calculated in the plenum. 

When the program is executed an estimation of the initial 
operating conditions must be supplied by the user. Besides 
velocity and thrust, initial conditions for plenum pressure, 
displacement and pitch angle must be supplied. From the 



displacement, pitch angle and velocity for calm water 
conditions, a plenum volume is computed. Using the leakage 
areas supplied, exit flow rates are computed for the bow 
and stern seals. The fan subroutine is then called using 
the pressure in the plenum, bow and stern seals as primary 
arguments. Using these pressures supply flow rates are 
established for. the fans that supply air to the plenum, 
and both seals. As the seals discharge into the main plenum 
the sum of these fan flow rates establish the total influx 
rate. The net flow being the difference between the influx 
and efflux rates establishes the new plenum pressure for- 
the next calculation interval based upon the adiabatic 
perfect gas law. This new pressure creates a new plenum 
lifting force which is then balanced by the defining differ- 
ential equations of motion resulting in a new operating 
condition of craft attitude and displacement. 

In studying craft motion in calm water the leakage area 
is of little or no consequence since no additional leakage 
areas are normally developed beneath the sidewalls or seals, 
and the rate of change of plenum volume is relatively slow. 
On the other hand, when the program is executed in simulated 
sea states where gap areas are generated beneath the seals 
and sidewalls and plenum volume does change quite rapidly 
the leakage areas and fan maps become of prime importance. 

It must be emphasized though that it is not the leakage 
area itself that affects craft performance. It is rather 
the relative size of the fixed leakage area with respect to 
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the variable leakage areas generated around the periphery 
of the craft by wave motion and craft attitude and position 
which is important. 

B. FAN MAP ANALYSIS AND DERIVATION 

Isolation of the factors which caused the large 
discrepancy between observed steady state gap area and 
calculated area required an examination of both the efflux 
and Influx equations for possible errors. 

The equation used to establish leakage flow rates under 
the seals and the sidewalls are the same. This equation 
is given below 



Q = C^Aj^(2Ap/p)^/^ 

It’ is the general flow equation for a sharp edged orifice 
in turbulent flow conditions. As the orifice coefficient, 

C^, is the only variable in the equation that could be 
dependent on craft geometry and its range was limited between 
unity and six tenths; the efflux equations were eliminated 
as the source of the error. The only area left for possible 
error was the fan maps themselves. 

The fan maps in question are arrays which give fan flow 
rates in cubic feet per second as a function of plenum or 
seal pressure. They are entered as data in "INCON" as 
characteristics of the specific craft being modeled. 
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It was found that the fan maps constructed by Leo and 
Boncal to represent the fans installed in the XR-3 were based 
upon incorrect information. Actual fan performance data for 
the fans currently Installed is presented*- in Figure 1, and 
physical characteristics in Figure 2. 

The fan performance map given by the manufacturer, Joy 
Manufacturing Company, is based upon flow that is not impeded 
by air ducts. It is therefore optimistic for any system 
that contains long ducts as these ducts will modify 
considerably the fan performance by creating ducting losses. 

The best way to obtain the actual fan performance maps 
in a duct system would be to actually measure the flow rates 
as the pressure in the receiving reservoir (in this case the 
plenum) . is varied. As this method was impractical in the 
case of the XR-3 another method of estimating the modified 
fan performance had to be devised. 

The method outlined below was utilized and although it 
probably does not give exact results the fan maps thus 
achieved should adequately represent the actual fan performance 
of the XR-3. 



^The fan data was provided by Associate Professor of 
Aerology D.M, Layton who is in charge of the actual XR-3 
test program at the Naval Postgraduate School. He has 
verified that the fans Installed are Joy Model AVR87-55D2058 
for which the data given in P’igures 1 and 2 are correct. 
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FIGURE 1 



Performance and Fan Map Curves for the Joy Ax i vane Fan 
(Model AVR 87-55D2058) 
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Pressure Inches W.G. — 6.5 Design Specification — Military 
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FIGURE 2 . . 

Manufacturer's Specifications for the Joy Axivane Fan 



The first step used In arriving at the new fan maps was 
to establish a realistic leakage area. Once this leakage 
area was established a steady state efflux rate was 
established. Because in steady state operating conditions 
the influx and efflux rates must be equal the flow rate 
supplied by the lift fans in total was established. It was 
at this point that a second assumption was made. This 
assumption was that there were greater losses in the ducting 
supplying air to the bow seal and main plenum chamber than 
there were in the ducting supplying the rear seal. This 
assumption was based upon the fact that the ducting to the 
bowseal and main plenum had to carry the air forward 
approximately sixteen feet to their discharge orifices 
whereas the stern seal ducting was only a few feet long with 
only one curve in it to impede flow. 

As a result the fraction of total flow assigned to the 
stern seal was greater on an individual fan basis than was 
the air supplied by the bov; seal or main plenum fans. 

Another assumption was that the same map could adequately 
represent both the bow seal supply fans and those supplying 
the main chamber. This last assumption was based on similar 
duct lengths and construction. 

Now knowing the individual fan flow rates and the actual 
plenum pressure a pressure drop in the duct for steady state 
conditions was established. This pressure drop was the 
difference betv;een the known plenum operating pressure 
(24.84 Ibs/ft^) and the indicated fan discharge pressure 
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for that flow rate. This pressure drop v;as then scaled 
for different flow rates based upon the fact that pressure 
losses vary with the square of the flow rates. In this 
manner new fan maps were computed for the bow seal and 
main plenum, and the stern sea. These maps are presented 
in Figures 3 and 

The last step was to go back to the leakage area and 
make the slight adjustments needed to give the required 

p 

plenum steady state pressure of 2^1. 8^1 Ibs/ft . The final 

2 

leakage areas established were .08 ft for the bow seal 
2 

and 3.79 ft . for the stern seal. 

In the original fan subroutine the horsepower generated 
by the fans is computed based upon the influx rate and the 
pressure differential between the atmosphere and the plenum. 
As the actual work done per unit time by the fans this 
figure can be misinterpreted as the actual power needed 
by the fans to generate printed work rate. This would only 
hold true if the fans operated with one hundred percent 
efficiency and there were no ducting losses. In actuality 
fans of this nature v;hlch supply relatively large volumes 
of air at lov; pressure rarely operate at over seventy percent 
efficiency and can achieve that only over an extremely 
narrow operating range. These physical characteristics of 
fans were v;ell presented by R.D. Moyer in 19^6 in Ref, . 

For the fans that are used in the XR-3 the Joy Manufac- 
turing Company has provided the fan's horsepower requirements 
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Bow Seal and Plenum Chamber Fan Maps^ adjusted for ducting losses 
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on the same performance graph as the fan maps. This power 
requirement is plotted for the entire operating range of 
the fauis and is fairly consistent at about 2 horsepower. 

This horsepower requirement does not change under the new 
fan maps since the fan map discharge pressures at the fan 
exhaust remains the same. This horsepower map was entered 
in the fan subroutine as a table of look up and actual power 
requirements were computed from it using flow rate as the 
entering argument. 

In the COLPIL subroutine (the routine that prints out 
all data and draws all graphs) the fan power label was 
changed to "fan power (ideal)" to avoid possible confusion. 
In the fan subroutine the print option was modified to give 
the actual power required as well as the ideal power. Also 
the fan efficiency is computed and is given as an output of 
this subroutine. 

C. SIMULATED PERFORMANCE COMPARISONS BETWEEN 

ORIGINAL AND NEW FAN MAPS 

As stated earlier it was believed that changing the 
leakage areas and the fan maps would have a significant 
effect only when the computer model was run under conditions 
where the plenum volume underwent comparatively rapid 
fluctuations and/or the leakage areas around the periphery 
of the hull varied. 

These assumptions v;ere confirmed by simulated runs under 
varying sea conditions. In all these comparison tests the 
two XR-3 models used v/ere identical except that the original 
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fan maps were replaced by the newly developed ones in one 
model. Both simulated craft were run using the original 
seals as developed in Ref. 3. Craft weight was 5813 lbs, 
center of gravity 9.57 feet from the stern, speed was held 
constant at 20 knots. All runs were made into ahead seas. 

The sea conditions chosen were Irregular sea state 1, 

Irregular sea state 2, regular seas 1 and regular seas 2. 

The latter two were so named because the amplitudes and 
frequencies chosen for these regular seas runs were the major 
component of Irregular sea state 1 and irregular sea state 2 
respectively. The actual amplitudes and frequencies for 
these runs are tabulated on page 27 . 

In the simulation runs Into regular seas 1 dramatic 
differences were noted In the fluctuation of plenum pressure 
and heave accelerations. With the original fan maps, 

2 

variations In plenum pressure were limited to about 2.2 Ibs/ft 
whereas with the new maps the variations were over 6.9 Ibs/ft . 
This represents an Increase In fluctuations of over 315 
percent. Because the plenum normally provides over seventy 
percent of the supporting force for the craft similarly 
large variations were noted in heave acceleration both at 
the bov/ and at the center of gravity. Pitch motion, though, 
did not change much being less than .1 degree In both cases. 

Results for the runs in other sea conditions are 
tabulated on page 26 • The differences noted for the runs 
In Irregular sea state 2 and regular seas 2 v;ere not as great. 
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This has been attributed to the craft's ability to conform 
its attitude to the wave slope in the case of the longer 
wave lengths associated with greater sea states. This ability 
to conform then reduces the high frequency volumetric changes 
in the plenum chamber, thereby reducing pressure fluctuations. 
Figures 5 through 12 are the results of running both models 
in the simulated sea state one. By visually comparing the 
corresponding variances of the two models it is readily 
apparent that the fan maps do significantly change the 
performance characteristics of the model. By examining the 
table on page 28 it is also apparent that one must be careful 
in predicting behavior as the variation in measured variables 
does not vary directly with the sea conditions. The new fan 
maps have in effect defined a similar but different craft. 
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Run 



Run 



Run 



Run 



TABLE I. SEA STATE WAVE COMPONENTS 



1 Regular Seas 1 Heading 180 



Wave 


Components 


Amplitude (ft) 


Frequency 


(rad/ sec ) 




1 


.0598 


4.1526 


! Regular Seas 2 


Heading I80 






Wave 


Components 


Amplitude (ft) 


Frequency 


(rad/sec) 



1 .2721 1.8447 



! Irregular Sea 

Wave Components 


State 1 Heading I80 

Amplitude (ft) 


Frequency (rad/sec) 


1 


.0029 


2.553 


2 


.0324 


3.256 


3 


.0599 


4.1526 


4 


.0558 


5.2961 


5 


.0401 


6.7545 


6 


.0262 


8.6145 


7 


.0165 


10.9867 


8 


.0102 


14.0122 



1 Irregular Sea 

Wave Components 


State 2 Heading I80 

Amplitude (ft) 


Frequency (rad/sec) 


1 


.0168 


1.1786 


2 


.1464 


1.4745 


3 


.2721 


1.8447 


4 


.2668 


2.3077 


5 


.2052 


2.8871 


6 


.1410 


3.6118 


7 


.0928 


4.5185 


8 


.0600 


5.6528 
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TABLE II SIMULATED XR-3 PERFORMANCE WITH DIFFERENT 

FAN MAPS 



Comparison of XR-3 Loads and Motion Computer Program 
performance with original for maps with XR-3 with new 
fan maps are given below: 



RUN 


QUANTITY 


XR-3 


XR-3 






(original) 


(new) 


1. Regular 


Plenum Pressure 


24.8±1.09 PSF 


24.8±3.i<5 PSF 


Seas 1 


Pitch Angle 


3±.05 deg. 


2.7±.05 deg. 




Draft 


6.25±.05 In. 


5.9±.l in. 




C.G. Acceleration 


±.0^ g 


±.15 g 



2. Irregular 


Plenum Pressure 


24.8+2.0 PSF 


24.8±15 PSF 


Sea State 1 


Pitch Angle 


3.2±.l deg. 


2.5±.2 deg. 




Draft 


6.5±.05 in. 


5.75±.l in. 




C.G. Acceleration 


±.07 g 


±.4 g 



3. Regular 


Plenum Pressure 


23.75±5.75 PSF 


24.8±7 PSF 


Seas 2 


Pitch Angle 


5.1±4.75 deg. 


4.2±3.8 deg 




Draft 


7.25±2.3 in. 


6.0±2.9 in'. 




C.G. Acceleration 


±.23 g 


±.27 g 



4 . Irregular 


Plenum Pressure 


5 to 33 PSF 


-5 to 40 PSF 


Sea State 2 


Pitch Angle 


.8+2.2 deg. 


.8±2.5 deg. 




Draft 


8±6 in. 


10+7 In. 




C.G. Acceleration 


±.46 g 


±•.7 g 




o 



g 



3C2 






T irae 



QCG 

(seconds) 



::o 



FIGURE 5 

Computer Simulation Plenum Pressure verses Time^ Speed 20 knots 
Original Fan Maps> Sea State 1 
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FIGURE 6 

Computer Simulation Pitch Angle verses Time^ Speed 20 knots 
Original Fan Fiaps^ Sea State 1 . 
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FIGURE 7 

Computer Simulation C.G. Draft verses Time^ Speed 20 knots 
Original Fan Haps^ Sea State 1 
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Computer Simulation C.G. acceleration verses Time, Speed 20 knots 
Original Fan Maps, Sea State 1 
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FIGURE 9 

Computer Simulation Plenum Pressure verses Time^ Speed 20 knots 
New Fan Haps^ Sea State 1 
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FIGURE 10 

Computer Simulation Pitch Angle verses Time^ Speed 20 knots 
New Fan Maps> Sea State 1 
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FIGURE 11 

Computer Simulation C.G. Draft verses Time., Speed 20 knots 
New Fan Haps., Sea State 1 
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FIGURE 12 

Computer Simulation C.G. Acceleration verses Time^ Speed 20 knots 
New Fan Maps^ Sea State 1 
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Ill, SEAL MODEL REFINEMENTS 



A. INTRODUCTION AND BASIC MODELING TECHNIQUES 

Through tests conducted aboard the XR-3 ty the NFS test 
group headed by Associated Professor D.M. Layton of the 
Aeronautics Department It has been demonstrated that the 
forces generated by the seals have an appreciable effect on 
the overall performance of the test craft. By varying the 
lower travel limit of the bow seal and contouring Its shape 
the steady state pitch angle can be changed as well as the 
thrust versus velocity characteristics. Although not as 
pronounced, changing the support cable length of the stern 
seal also affected craft trim conditions. 

Because of the demonstrated steady stats pitch sensitivity 
of the XR-3 to seal position limits as set by the support 
cables It was decided that more detailed modeling ‘of the 
seals v/as needed before the overall response of the L&M 
computer model would reflect the behavior of the XR-3 under 
varying conditions of support cable length. 

The original XR-3 seal modeling as presented In Ref. 3 
was patterned after the rear seal in the L&K computer model 
of the Bell (lOOB) 100 ton CAB. This seal Is of the "soft" 
pneumatic type and Its representation is well covered by 
Kaplan, Bent son and Sargent In Ref. 2. 

The pneumatic "soft" seal model is based on the assumption 
that this type of seal is unstiffened and non-rigid. The 
mass of the seal v;as considered so insignificantly small as 
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compared to the forces generated by the seal when It Is In 
contact with the water as to have negligible effect on overall 
craft dynamics. The lifting forces generated by the seals 
were calculated by multiplying the wetted contact area of 
the seals by the pressure differential between the seal and 
the atmosphere for the bow seal and the plenum chamber and 
seal in the case of the stern seal. 

Drag forces on the seals were considered to be due only 
to skin friction between the water and the seal. The equation 
used to compute these forces was 

P = - C,.. A 

X T 2 di seal 



where p is the density of water, u the velocity of the craft 

and A T is the wetted surface of the seal. C,„ is the 
seal df 

drag coefficient and is calculated using the Schultz-Griinow 
formula 



Cdf = .^27/ 

where Re is the Reynolds number computed based on the average 
wetted length of the seal and the Reynolds number is 

calculated by : 

D _ u X (wetted length) 
ne F 

1.28 X 10“^ 

v;here 1.28x10“^ ft^/sec is the kinematic viscosity of the 



water. 



The actual construction of the seals is well described 
in Ref. 3. Each seal consists of two pneumatic lobes 
separated by a common membrane. The membrane has limber 
holes in it to equalize the pressure between the lobes. The 
lov/er surface of the bottom lobe is shaped by twelve steel 
battens that act as stiffeners. These battens are spaced 
uniformly across the bottom surface at right angles to the 
hinge line. It is this surface that comes in contact with 
the water, and will hereafter be referred to as the "leading 
face" of the seal. 

The forward edge of the leading faces is hinged near. the 
top of the plenum chamber allowing the seal to pivot about 
this line when waves are encountered. 

The pressure in the seals is maintained by axial fans. 

The air is delivered by ducts that discharge into the top 
of the upper lobe. The seals in turn vent to the main 
plenum by means of ducting. Back pressure in the vent 
ducting can be controlled by butterfly valves that can be 
manually adjusted. Normally the seals are operated so that 
they maintain about 1-2 PSF more pressure than the plenum. 

Dovmward travel limits of the seals are controlled by 
adjustable support cables attached to the rear edge of each 
steel batten at the trailing edge of the seal's leading face. 
In addition a second support cable is attached to each 
stiffener approximately midv/ay between the hinge line and 
the trailing edge. Figure 13 shov;s a cutav/ay of the seals, 
only two cables are drawn in for clarity. 
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FIGURE 13 

Cross Sectional View of Pneumatic Seal as installed on the XR-3 
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Both the bow and the stern seal extend the width of 
the plenum chamber and are identical in construction. The 
leading faces are approximately i»9 Inches long from hinge . 
line to trailing edge. The perpendicular distance from 
sidewall keels to the hinge line is 22.5 inches. 

As stated before the original modeling of the seals did 
not take into account any of the seal’s ability to flex or 
pivot on its hinge. In addition no attempt was made to 
model the support cables and their ability to affect craft 
attitude. In the new modeling the support cables have been 
Incorporated. In the stern seal modeling only the rear 
support cable is modeled since this seal rarely comes in 
contact with the water. . In the bow seal subroutine both the 
rear support cables and the center support cables have been 
included since this seal does contact the water. Before any 
attempt at revising the seal subroutines could be undertaken 
it was necessary first to study the action of the seals while 
the craft was actually underway. 

The semi-rigid construction of the seal dictated a highly 
complex distributed parameter model of the seals be developed 
if the seals were to be fully modeled without any approxima- 
tions. In addition to the fact that practically all of the 
physical characteristic parameters of the seal such as spring 
constants for the steel stiffeners and fabric flexibility 
v/ere unknovm such a model would require an inordinate amount 
of calculation time in the computer. Since extra computation 
time v?as to be avoided simplifying assumptions had to be made. 



To be accurate these assumptions or approximations had to 
be based on actual observations of the seal or on previously 
made approximations whose effect has already been investigated. 

In Ref. 2 the approximation of eliminating seal dynamics 
and the setting of seal pressure at a constant level above 
plenum pressure has been thoroughly investigated and was 
incorporated in the original model of soft pneumatic seals. 
These approximations have been continued in the present seal 
modeling. 

The other modeling approximations made were based upon 
actual observations made of the seals while the XR-3 was in 
motion. These observations are described below: 

1. For small v;ater perturbations (approximately 3 Inches) 
the bcv; seal flexed and appeared to pivot slightly about the 
center support cable causing the angle at the hinge line 

to slightly Increase, but almost Insignificantly. 

2. For longer perturbations the seal both flexed as 
described above and then began to pivot on the hinge. After 
the seal rose about six Inches no additional flexing took 
place . 

3. The bovj seal maintained contact with the water at all 
times and never gaped because of an inability to follow wave 
motion . 

The bow seals penetration of the water was practically 
independent of water motion in relation to the craft, but 
appeared to be inversely related to the craft's velocity. 
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When the XR-3 was at rest penetration was about five to six 
inches and when traveling at 20 knots appeared to be 
approximately 1 inch. 

5. The stern seal rode out of the water while the craft 
operated above "hump" speed. The height of the air gap thus 
created appeared to remain almost constant at four to five 
inches. Water rarely contacted the seal. 

Using the above described observations and the preceding 
simplifying approximations concerning seal dynamics and 
pressures both the bow and stern seal were modified to 
reflect the action of the actual seals. 

The approximation given in Ref. 2 of equating the wetted 
length of the seal with the water height above the sidewall 
keel at the seal's trailing edge location was eliminated. 

In its stead a two-dimensional "table of look up" was 
developed which gave wetted length as a function of the 
immersion depth of the seal's trailing edge in the water and 
the seal's trailing edge deflection height above the keel at 
the trailing edge location. In other words if the seal was 
deflected upwards because of wave motion or pitching motion 
of the craft the deflection height of the trailing edge 
would increase and the seal's leading face would become more 
parallel to the water's surface thereby increasing the wetted 



area of the seal. 



The "table of lookup" was constructed using the 
information obtained from a graphical deflection model 
similar to the one presented in Figure 14. In this model 
the leading face of the seal is represented in varying 
amounts of deflection height from full down (point A) to 
full up against the seal framework (point B). Full down 
is zero deflection and is where the trailing edge of the 
seal is even with the bottom of the sidewall keel, full- up 
is maximum deflection. This deflection height was divided 
into segments of one inch for ease of computation. Points 
were marked upon each leading face line which represented 
the surface of the water on the leading face under varying 
depths of trailing edge immersion. The distance along the 
leading face to the trailing edge from these points then 
represented the length of the leading face that was immersed 
in the water or the "wetted length". The points for sim.ilar 
immersion depths v;ere connected for the varying seal deflec- 
tion heights and these construction lines are labeled with 
their corresponding immersion depths in Inches along the 
zero deflection height leading face. The construction line 
corresponding to an immersion depth of zero inches is, of 
course, congruent with the trailing edge of the seal since 
if the seal did not penetrate the water to at least some 
degree there could be no "wetted length" unless the leading 
face was exactly parallel to the water’s surface, which 
never occurs. These construction lines v;ere dravm for 
Immersion depths of the trailing edge of 0 to 5 inches. 
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FIGURE 14 

Bow AND Stern Seal Deflection Model ( not to scale 
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The actual immersion depth of the trailing edge is 
calculated based upon the observation that at rest the 
immersion depth was approximately 5.5 Inches and while 
traveling at 20 knots was estimated at a little less than 
one inch. Since the hydrodynamic force on the seal is 
proportional to the square of the craft's velocity the 
Immersion depth was assumed to be approximately inversely 
proportional to some factor of the square of the speed. 

Using 5.5 Inches as the known immersion depth at zero 
zpeed and .75 inches at 20 knots, the following empirical 
formula was derived which gives reasonable results 

D(in.) = 5.57(1 + 

where D(in.) is the depth in Inches and u is the velocity 
in feet per second. 

Using the calculated Immersion depth and the" seals 
calculated deflection height the table is entered to give 
the wetted length. This table of look up is presented 
graphically in Figure 15. The non-linear nature of the 
curves is readily apparent. An Important aspect of this 
table apparent in Figure 15 is that within about the first 
ten inches of seal deflection the v;etted length of the seal 
remains fairly constant, especially if the craft is operating 
at speeds v/here the immersion depth is approximately one inch 
This figure also shov/s that the maximum v/etted length is 
equal to the length of the leading face (just over 48 inches) 
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Wetted Length (inchetj) 

FIGURE 15 

Seal Trailing Edge Height above Keel verses Seal Wetted Length 



As the seal deflects and becomes more horizontal the wetted 
surface area becomes greater given a constant Immersion 
depth. 

This same "table of look up" Is provided In the stern 
seal modeling but because the seal does not ride In the water 
an assumed "depth" had to be used as an entering argument. 

The assumption was made that since the stern seal does 
not ride on the water Its lifting force has to be the result 
of some form of air pressure. This force could only be 
generated by the impingement of air particles on the seal's 
leading face as they rushed to escape underneath the seal. 

At this point an assumption was made that the effective 
wetted area occurred where the kinetic energy in the air mass 
was approximately one quarter that of the air escaping under- 
neath the seal. As the air in the plenum is at such a 
pressure that it can be considered incompressible this would 
Occur where the cross-sectional area was approximately two 
times the cross-sectional area beneath the seal. As the 
width of the plenum chamber is constant and the gap height 
is approximately 4,5 inches at all times the point on the 
seal where the cross-sectional area is two times the gap 
area is the point along the leading face where the height 
from the seal to the water is two times the gap height or 
about 9 Inches. Subtracting the gap height gives an 
effective "phantom" Immersion depth of 4.5 inches; this was 
rounded to five inches. 



Lift forces generated by the stern 'seal were considered 
to be the product of the seal area found by using the above 
phantom immersion depth and the seal deflection height, and 
the pressure differential between the stern seal and plenum 
chamber. 

The above lifting force calculation has only minimal 
effect on simulated performance of the XR-3 since the 
pressure between the seal and the plenum chamber is normally 
only one to two PSP. The drag forces on the stern seal are 
calculated in the same manner as they are in the bow seal 
subroutine but are subsequently set to zero since the seal 
rarely contacts the water. 

B. BOW SEAL MODELING 

Contained herein is a detailed analysis of the actual 
bow seal subroutine. It is coded in standard FORTRAN. By 
reviewing the actual subroutine the actual approximations 
used in the routine may be understood. A full listing of 
the BOVJSL subroutine is contained in Appendix A along with 
the required entries in the data input statements. 

(1) DO 1 J=l,ll 

(2) GAP(J) = 0.0 

(3) ELSKI(J) = 0.0 

{h) WETLEN(J) = 0.0 

(5) ELSKID(J) = 0.0 

(6) 1 CONTINUE 

The above statements (1 through 6) zero all variable arrays 
that will be used every time BOVJSL is entered. GAP(J) is 
the height of the air gap at each station along the bow seal. 



ELSKI(J) is the water height above the keel at the trailing 
edge of the seal. WETLEN(J) is the wetted surface length 
at each station. ELSKID(J) is the deflection height of the 
seal trailing edge measured from the botton of the keel of 
the craft. 

(7) ALBS =0.0 

(8) FX = 0.0 

(9) PY = 0.0 

(10) PZ = 0.0 

(11) PK = 0.0 

(12) PM = 0.0 

(13) PN = 0.0 

The preceding seven statements initialize the single valued 
variables that are computed each time BOWSL is entered. 

ALBS is leakage area of the bow seal. PX is the forces 
parallel to the craft's axis of symmetry (centerline). 

FY is the horizontal forces perpendicular to the FX forces. 

FZ is the vertical forces. FK, FM, FN are the moments about 
the Z axis, Y axis and X axis respectively. 

(14) DELPBG = PBS ~ PB 

(15) IF (DELPBG.lt. 0.0 )DELPBG=0.0 

Equation l4 computes the pressure differential between 
the plenum chamber (PB) and the seal (PBS). This differential 
is actually fixed in INCON as an entering argument. Equation 
15 assures that the seal differential pressure is never 
negative. This is reasonable since the seals vent to the 
chamber. 
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(16) ARGO = ELMAXB/CORLEN 

(17) ANGO = ARSIN(ARGO) 

(18) XI = XBS+ZBS*THETA-CORLEN*COS(ANGO) 

The three equations above (I6-I8) establish the moment 
arm in the X direction from the craft center of gravity 
to the trailing edge of the bow seal. ELMAXB is the rear 
cable length. CORLEN is the chord length of the bow seal. 
XBS is the distance from the center of gravity to the bow. 
ZBS is the vertical distance of the bow hinge line below 
the center of gravity. Theta is the pitch angle of the 
craft . 



(19) Z1 = -Z-ABS+XBS5t'THETA-ELMAXB*C0S (THETA) 

Equation (19) establishes the moment arm in the -Z 
direction to the water's surface through which the drag 
forces act to create a moment about the Y axis. 

(20) DPTHFT = (5. 5/(1. +(U/25. ))**2. 0/12.0 

Equation (20) calculates the immersion depth of the 
seal based upon craft velocity, u. 

(21) IP(CENCAB.GT.l. 1875 )CENCAB=1. 1875 

Logical statement (21) sets the maximum length of the 
center support cable to 3 inches more than the maximum rear 
cable length, ELMAXB. 
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( 22 ) 



DO 3 K=1,N 



3 CONTINUE 

The above loop calculates for each seal partition point 
the water height above the keel ELSKl(K), the wetted length, 
WETLEN(K), and the seal deflection height ELSKID(K) in 
preparation to calculating the forces and moments. 

(23) DPFT(K) = DPTHFT 

Equation ( 23 ) initializes the immersion depth of the 
seal . 



(24) ELSK1(K)=(ETA(3,K)-DTABX(K)*(XX(3,K)-X1)-Z1) 

+YY ( 3 , K ) *PHI+XLXPWV*V/ATSLP 

Equation (24) establishes the water height above (or 
below) the keel at the bow. ETA(3,K) is the height of the 
water at the center of gravity above mean. DTABX(K) is 
the correction factor for moving this height away from the 
C.G. (Center of Gravity). (XX(3,K)-X1) is the distance the 
correction factor has to act on. Z1 is the height of the 
keel in relation to the C.G. YY(3,K) is the displacement 
of the partition from the craft center line. PHI is the 
roll angle. XLXPNV*WATSLP is the correction factor for the 
water slope caused by the plenum pressure acting for a 
finite period. 
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(25) 



IF(ELSKKK) .GT.HINGHT) ELSKI (K) =HINGHT 



Equation (25) limits the water height to the hinge 
height of the bow seal. 

(26) IF((HINGHT-ELSKI(K)+DPPT(K)) .GE.ELMAXB) 
DPFT(K)=ELMAXB-HINGHT+ELSKI(K) 

Equation (26) adjusts the immersion depth of the seal 
so that the seal can go no lower than that permitted by the 
rear cable ELMAXB. 

(27) IP(DPFT.LT.O.O) DPFT=0.0 

DPFT is used in computing DPIN which is an entering 
argument into the wet length table BWSL. Therefore 
equation ( 27 ) limits this argument to a minimum value of 
zero. 

(28) ELSKID(K)=(ELSKI(K)-DPFT(K) )*12 .0 

Equation (28) calculates the deflection height of the 

seal (ELSKID(K) ) in Inches based upon water height, ELSKI(K) 
and the immersion depth DPPT(K). 

(29) IP( (HINGHT-ELSKID(K)/12.0) .GE.ELMAXB) 

ELSKID(K)=(HINGHT -ELMAXB) ^12.0 

Equation (29) limits the downward travel of the bow seal 
to ELMAXB the rear cable length and adjusts the deflection 
height ELSKID(K) to reflect this lower limit. 
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(30) • DPIN=DPFT(K)*12.0 

(31) MM=DPIN 

(32) MM1=MM+1 

(33) MM2=MM1+1 

(3^) DINC=DPIN-MM 

The equations above, (30) through (3^)j establish 
parameters for entering the bow seal wetted length array, 
BOWSL. DPIN is the immersion depth in inches, MMl is the 
lower limit table calling integer, MM2 is the upper table 
calling integer and DING is the incremental difference. 

( 35 ) gap ( K ) =-ELSKI (K ) + (HINGHT-ELMAXB ) 

Equation (35) computes the air gap height based on 
water height ELSKI(K) which is positive upwards and the 
limit of the seal travel (HINGHT-ELMAXB). 

( 36 ) IF(GAP.LT.O.O) GAP= 0.0 

Equation ( 36 ) limits the gap to positive values. 

(37) IF(ELSKID(K) .GE.0.0) GO TO 2 

( 38 ) WETLEN(K)=ELSKI(K) 

(39) GO TO 3 

Equations (37) through (39) set the wet length equal to 
the water height if it is negative. This computation is 
used in case one seal partition point has a negative value 
while the others on either side are positive. 
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(110) MM3=ELSKID(K) 

(111) MMi}=MM3+l 

(i}2) MM5=MMil+l 

(il3) DLING=ELSKID(K)-MM3 



Equations (i|0) through (il3) establish the final two 
entering arguments for retrieving the four points needed 
from BOWSL to compute the wetted length. ELSKID(K) is the 
deflection height of the seal. MMil is the lower entering 
integer based upon the deflection height , MM5 is the upper 
Integer. DLINC is the incremental difference used for linear 
interpolation. 



(ilil) BWSLl=BWSL(MMl,MMll) 

(115) BWSL2=BWSL(MM1,MM5) 

(116) BWSL3=BWSL(MM2,MMil) 

(117) BWSLl}=BWSL(MM2,MM5) 

( ll 8 ) BWSLA1= ( BWSL2-BWSL1 ) *DLINC+BWSL1 

(Il9) BWSLA2=(BWSLil-BWSL3)*DLINC+BWSL3 

( 50 ) WETLEN (K ) = ( ( BWSLA2-BWSLA1 ) *DINC+BWSLA1 ) /12 . 0 



Equations (il^l) through (50) identify the four points of 
the BWSL array which bracket the actual value of wetted 
length. Linear interpolation is used twice to find the 
average based on deflection height and once to find the 
final wetted length. The array BWSL was calculated in inches 
thug, to convert to feet the final average is divided by 
12 . 0 . 



(51) DO 8 J=1,N 



8 CONTINUE 

This loop makes all the force calculations based upon 
wetted length and seal position calculated previously. In 
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addition it calculates the leakage area based upon the steady 
state leakage area and the variable gap heights. 

(52) WETLAV=(WETLEN(J+1)+WETLEN(J) )/2.0 

(53) IF(WETLAV.LE. .001)00 TO 6 

(5^) DPFTAV=(DPFT(J+1)+DPFT(J) )/2.0 

(55) ELSKIA=(ELSKI(J+l)+ELSKI(J))/2.0 

(56) ELSKDA=(ELSKID(J+1)+ELSKID(J) )/2.0 

Equations (52) and (5^) through (56) average the 
wetted length, immersion depth, water height and deflection 
height respectively, for each seal segment based on the 
values calculated for the stations on either side. 

Equation (53) routes the program around the remaining force 
calculations if the wetted length of the seal is zero or 
insignificant for calculation. 

(57) SEALHT=HINGHT-ELSKDA 

Equation (57) is the calculation of the distance betv;een 
the top of the seal's supporting frame and the bottom of 
the trailing edge of the seal. 

(58) . DIFF=2.0*CENCAB-(SEALHT+0.5) 

(59) IF(DIFF.GT.0. 5) DIFF=0.5 

The DIFF equation (58) calculates the spring distance 
of the steel battens based upon the center support cable 
length CENCAB and the position of the seal, SEALHT. 

Equation (59) limits the spring movement to .5 feet or 
six Inches, the point v;here rigid body analysis comes into 
full effect. 
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(60) ARMlB(J)=Xl+WETLAV/2.0 

( 6 1 ) ARM2B ( J ) = ZS-ELSKI A+DPFTAV/2 . 0 



Equations (60) and (6l) calculate the moment arms in 
the X and Z direction respectively based upon the length to 
the trailing edge of the bow seal, XI, and half the wetted 
length, WETLAV/2.0, for ARMIB. ARM2B is based upon the 
height of the center of gravity ZS , the water height, 

ELSKIA and half the Immersion depth of the seal, DPFTAV/2.0. 



(62) IF(DIFP.GE.0.25) GO TO 

( 63 ) DFBS ( J ) =-DELP ^•^DEL YB S * WETLAV 

(64) GO TO 5 



Equations (62) through (64) calculate the vertical force 
acting on the segment if the deflection difference is less 
than .25 feet or three inches. If the deflection difference 
is greater rhan .25 or 3 inches the vertical force equation 
has a spring constant term in it acting upon the moment 
summations about the hinge line. 



( 65 ) 4 FORLEN=XBF-WETLAV 

(66) IF(PORLEN.EQ.OO) GO TO 5 

( 67 ) ARGW= (HINGHT-ELSKI A ) /FORLEN 

(68) IF(ARGW.GT. 1.0)ARGW=1. 0 

( 69 ) ANGV;=ARSIN(ARGW) 

( 70 ) porcos=cosungv/) 

(71) GO TO 6 

( 72 ) 5 FORCOS-0.0 



Equations ( 65 ) through (72) calculate the length of the 
seal not in contact with the v;ater, FOKLEN, and the angle 
which this length makes with the water ANGW. From the angle 
the cosine of ANGV/ is, calculated, FORCOS, for future use 
in calculating the vertical moments about the hinge line. 



57 



(73) DFBS(J)=-DELP*DELYBS*WETLAV-DELP*FORLEN*DELYBS 

*FORCOS * ( FORLEN *0.5 *FORCOS ) / 
F0RLEN*F0RC0S+WETLAV/2 . 0 ) ) * 
((DIFF-0.25)*4.0) 

The above equation calculates the vertical force acting 
upon the bow seal by summing vertical moments about the 
hinge line and equating it to the moment created by one 
vertical force vector acting on the midpoint of the wetted 
surface. The term " (DIFF-0 . 25) . 0" introduces a spring 
distance of three inches through which the vertical force 
transitions from being only based upon the product of the 
wetted area and the pressure differential to one based upon 
equalizing moments about the hinge line. After the seal 
bends this three inches it is assumed to be completely rigid 
and therefore if the seal is not undergoing accelerations, 
the moments must balance about the hinge line. 

(7^0 ARG=0.5*RHO*U*U*WETLAV*DELYBS 

(75) reski=u*wetlav/enu 

(76) CDTSKI=0. i}2 7/ (ALOG10( RESKI )-0. 407)^*2.64 

(77) TSKIB(J)=-ARG*CDTSKI 

The above four equations calculate the drag forces of 
the seal based upon the previously discussed Schult z-Grunow 
formula. RESKI is the Reynolds number ,CDTSKI is the drag 
coefficient and TSKIB(J) is the drag force on the segment. 



7 CONTINUE 

(78) FX=FX+TSK1B(J) 

(79) FZ=FZ+DFBS(J) 

(80) FK=FK+DFBS(J)*YAVGB(J) 

( 81 ) FM=FM-DFBS ( J ) *ARM1B ( J ) +TSKB ( J ) * ARM2B ( J ) 

(82) FN=FN-TSK1B(J)*YAVGB(J) 

( 83 ) ALBS=ALBS+GAP(J)+GAP(J+l))*DELYBS/2.0 

Equations ( 78 ) through ( 83 ) calculate the forces and 
moments for each segment and sums them to give total forces 
and moments. YAVOB(J) is the average distance off centerline 
of the individual segment. 

(8ii) 8 CONTINUE 

( 85 ) ALBS=ALBS+BLEAK 

(86) SQFAC=SQRT(2. *ABS(PBAR)/RHOINF) 

( 87 ) QL=CFBS*ALBS*SQFAC*SIGN( 1 . ,PBAR) 

The last equations calculate the flow rate of air based 
upon the leakage area ALBS and the pressure differential 
between the plenum chamber and the atmosphere. BLEAK is 
the seal leakage entered in INCON to account for slight 
leakage vjhere the seals do not quite maintain contact with 
the sidewall. 

The remaining statements of subroutine BOWSL give the 
following values as printed output if the bov; seal print 
switch is on, GAP, WETLEN, FX, FY, FZ , FK, FM, and FN. 

C. STERN SEAL MODELING 

The stern seal subroutine is covered in the same manner 
as the bow seal subroutine. The complete listing of the 
stern seal subroutine is contained in Appendix B along v;ith 
instructions on the proper entries to be used on the data 
cards. 
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(1) DO 1 J=l,ll 

(2) GAP(J)=0.0 

(3) ELSKI(J)=0.0 

(4) ELSKIL(J)=0.0 

(5) AIRLEN(J)=0.0 

(6) 1 CONTINUE 



The above loop and equations zero all variable arrays 
that will change each time subroutine STNSL is entered as 
was the case in the bow seal subroutine. 



(7) EFDEP = 6.0 

Equation (7) establishes the effective depth of the seal 
at five inches. This figure is used as an entering argument 
in the stern seal contour mapping array CTNSL. 



(8) 


ALSS=0.0 


(9) 


FX=0. 0 


(10) 


FZ=0.0 


(11) 


FK=0.0 


(12) 


FM=0. 0 


(13) 


FN=0.0 


(14) 


AGAP1=0.0 


(15) 


AGAP2=0. 0 


( 16 ) 


AGAPA1=0.0 



The nine equations above zero all single valued variables 
that are computed each time STNSL is entered. 



(17) DELP=PSS-PB 

(18) IF(DELP.LT.O.O) DELP=0.0 

The two statements above fix the pressure differential 
betv/een the plenum (PB) and the stern seal (PSS). Equation 
(17) is first used to arrive at the difference between 
these tv.'o quantities and the logical IF statement is used 
to insure all values are non-negative. As in the bow seal 



• 60 



subroutine this Is logical since the air supplied to the 
stem seal vents Into the plenum chamber. 

(19) PBAR=PB-PINP 

This equation establishes the pressure differential 
between the plenum (PB) and the atmosphere. It Is used In 
computing leakage flow rates. 

(20) SINDIF=SINTH-COSTH*THETA 

(21) COSDIP=COSTH+SINTH*THETA 

The above equations (20) and (21) establish the angle 
of the XR-3 with respect to the horizontal and the vertical, 
respectively. 

(22) X1=XSS+ZSS*THETA-XLP«SINDIP 

(23) Z1=(-Z-ZSS+XSS*THETA-ELMAXS*C0S (THETA) ) 

The above two equations esteibllsh lengths In the X and Z 
directions which are used to calculate the water height above 
the keel at the stern seal location. 

(2^4) N=NSTA(i|) 

(25) DO 2 K=1,N 
2 CONTINUE 

The above loop Is used to establish the gap at each 
station and the phantom "wetted length" of the seal. 

NSTA(i|) is the number of stations chosen for the seal and 
can be up to eleven, thereby defining ten seal segments. 
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(26) ELSKI(K) = (ETA(4,K)-DETADX(K)*(XX(i|,K)-Xl) 

-Z1)+YY(4,K)»PHI 

(27) IF(ELSKKK) .GT.HINGHT) ELSKI(K)=HINGHT 

Equation (26) computes the height of the water above 
the keel at the stern seal location. ETA(i},K) Is a length 
correction factor that Is dependent on seal shape. Since 
In the case of the XR-3 the seal Is flat across the hinge 
line this term Is zero. DETADX(K) Is a wave correction 
factor that Is multiplied by the seal distance from the 
center of gravity (XX(4,K)-X1) to give the wave height at 
the stern seal location. The logical IF statement limits 
the water height to the top of the seal. 

(28) ELSKIL(K)=ELSKI(K)+GPS 

(29) IF(ELSKIL(K) .GT.HINGHT) ELSKIL(K)=HINGHT 

(30) IF(ELSKIL(K) .LT. (HINGHT-ELMAXS) ) 

ELSKIL ( K ) =HINGHT-ELMAXS 

ELSKIL(K) Is the height of the trailing edge of the 
stern seal above the keel. It is measured in feet. GPS is 
the air gap height that was observed to be constant under 
normal running conditions. It Is .4644 feet. The logical 
IP statement (29) establishes the upper limit of_seal 
deflection where It Is stopped by the seal supporting 
framev/ork. The second logical IF statement establishes the 
lower limit of seal deflection height. 

(31) GAP(K)=-ELSKI(K)+ (HINGHT-ELMAXS) 

(32) IP(GAP(K).LT.0.0)GAP(K)=0.0 

The equation (31) and logical IF statement (32) establish 
the gap height below the lov/er stop limit set by the support 
cable length ELMAXS. 
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(33) MM1=ELSKIL(K)*12.0 

(3^) MM2=--MM1 

(35) MM3=MM2+1 

(36) DLINC=ELSKIL(K)*12.0-MM1 

(37) STNSL1=CTNSL(MM,MM2) 

(38) STNSL2=CTNSL(MM,MT43) 

( 39 ) AIRLEN ( K) = ( ( STNSL2-STNSL1 ) *DLINC+STNSL1 ) /12 . 0 



Equations (33) through (35) establish the entering 
arguments of the seal contour array CTNSL In the same manner 
as was done In the bow seal subroutine. AIRLEN(K) Is the 
phantom "wetted length" and Is given the name AIRLEN 
because the seal Is supported by air. 



(40) N=NSTA(4)-1 

(41) DO 5 J=1,N 
5 CONTINUE 



The above loop accomplishes the actual force and 
moments of each seal segment and totals the forces and 
m.oments. It also computes the average gap height based upon 
the gap heights of the Individual stations. 

(42) ELSKIA=(ELSKI(J+1)+ELSKI(J) )/2.0 

(43) ELSKLA=(ELSKIL(J+l)+ELSKIL(J))/2.0 

(44) AIRLAV=(AIRLEN(J+l)+AIREN(J))/2.0 

The three above equations average the water height above 
the keel, ELSKIAjthe average seal deflection height ELSKLA^ 
and the average phantom wetted length AIRLAV. 

(45) AGAP=ELSKLA-ELSKIA 

(46) AGAP1=AGAP 

Equation (45) establishes the average height of the air 
gap at the seal segment. AGAPl Is a variable used In further 
calculations . 
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(47) IF(AGAP.LT.GPS)AGAP=GPS 

(48) IF(AGAP1.GT.GPS)AGAP1=GPS 

The above logical IF statements establish the lower 
limit for AGAP at .4644 feet at the maximum value of AGAPl 
at .4644 feet. 

(49) ARM1S=XX(4,J)+ELSKIA/2 

(50) ARM2S=ZS-ELSK1A 

The above equations (49) and (50) calculate the moment 
arms through which the vertical lifting force acts to create 
the major component of the FM moment about the Y axis; the 
second equation creates the moment arm with respect to the 
C.G. through which the drag forces act. 

(51) DFSS(J)=-DELP*DELYSS*AIRLAV/(GPS/AGAP)**2. 0 

Equation (51) calculates the vertical force acting upon 

the seal and with the term (GPS/AGAP) **2 . 0 allows for this 
force to decrease with the square of the ratio of the actual 
gap v;ldth AGAP and the nominal width GPS. This is based on 
the assumption that the air exiting beneath the stern seal 
supports the stern seal and that the lift forces generated 
by the exiting air are inversely proportional to their 
velocity. 



(52) ARG=.5RH0*U*U*AIRLAV*DELYSS 

(53) RESKI=U*AIRLAV/ENU 

(5^1) CDTSKI=. il27/(ALOG10( RESKI )-. 407) **2.6il 

(55) TSKIS(J)=-ARG*CDTSKI 

(56) TSKIS(J)=0.0 



Equations (52) through (55) calculate the theoretical 
drag forces based upon the average phantom wetted length 
AIRLAV as if the seal was actually in contact with the water. 
The Schultz Grunow formula is used as in the bow seal routine 
to calculate the drag forces based upon the Reynold's number 
which is based on AIRLAV the average "wetted” length. 

Equation (56) removes this drag calculation since the seal 
is presumed to ride on a cushion of air and is hence 
negligible . 



(57) FX=FX+TSKIS(J) 

( 58 ) FZ=FZ+DFSS(J) 

(59) FK=FK+DFSS(J)*YAVGS(J) 

( 60 ) FM=FM-DFSS ( J ) *ARM1S ( J ) +TSKIS( J ) *ARM2S ( J ) 

( 61 ) FN=FN-TSKIS(J)*YAVGS(J) 



The five equations above calculate the forces and 
moments of each segment and sums them for all the segments 
to give the total seal forces and moments. 



( 62 ) ALSS=ALSS+(GAP(J)+GAP(J+1) )*DELYSS/2.0 

( 63 ) AGAP2=AGAP2+AGAP1 

(64) AGAPA1=AGAP2/J 



The above equation (62) calculates the total gap area 
ALSS that is due to the seal reaching the downward limit 
of its travel and the gap thus increasing past its nominal 
value of .4644 feet. Equations ( 63 ) and (64) calculate the 
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average gap height including the nominal height and 

allows the nominal gap to close If the seal is slammed 
against its upper limit, the upper seal support framework. 

( 65 ) ALSS=ALSS+ALEAK* ( AGAPAl/GPS ) 

Equation (65) calculates the leakage area under the seal. 
ALSS is the leakage area due to the gap width increasing 
beyond the nominal .H6HH feet. The term ALEAK is the 
nominal leakage area of 3*79 square feet normally observed. 
Multiplying ALEAK by AGAPAl/GPS allows this nominal leakage 
area to decrease if the seal is limited in its upward travel. 

(66) SQFAC=SQRT(2. »ABS(PBAR/RHOINF) 

( 67 ) QL=CFSS*ALSS»SQFAC«SIGN(1. ,PBAR) 

Equations (66) and ( 67 ) calculate the flow rate based 
upon the leakage area and the pressure differential between 
the plenum and the atmosphere. The equation assumes orifice 
type leakage. 

The remaining statements give print outputs of the terms 
GAP, AIRLEM, FX, FY, FZ , FK, FM and FN if the stern seal 
print sv;itch is on. 
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IV. DATA ANALYSIS OP ACTUAL CRAFT TRANSIENT 
MOTION AND COMPUTER SIMULATION 



A. INTRODUCTION 

In order to verify the modeling changes made In the 
seals, data was collected on the XR-3 test craft for 
conditions of Induced transient motion In pitch and roll. 

To Induce roll action the test craft was sallied by 
shifting body weight laterally near the C.G. until maximum 
oscillations built up, then the forcing action was stopped. 
For pitch oscillations the chase boat was veered across 
the path of the XR-3 to produce a wave front with Its wake. 
By having the XR-3 pass over the wake It was Intended to 
Induce a transitory pitch oscillation. 

These attempts to produce a transient response cannot 
be vlev;ed as entirely successful. The principal reason for 
the limited success was that It was extremely difficult to 
Induce a large enough perturbation In either pitch or roll 
so that the normal background oscillations could be Ignored. 
Figures l6, 17 and l8 are graphs of the roll and pitch 
motions of the XR-3 when making Its twenty knot calibration. 
The magnitude and frequency of the steady state pitch 



^Figure 17 depicts the actual pitch angle of the test 
craft. Figure l8 Is the same data with the steady state 
average pitch angle removed. This was done on all data to 
facilitate the measurement of the pitch oscillations super- 
imposed on the steady state value of pitch angle. 
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FIGURE 16 

Roll Motion verses Time^ Velocity 20,^29 knots 
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FIGURE 17 

Pitch verses Time, Velocity 20.^29 knots 
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oscillation shown in Figure l8 is approximately .08- degrees 
and .295 Hz respectively. The steady state magnitude of 
the roll motion shown in Figure I 6 was approximately .075 
degrees and its frequency was approximately the same frequency 
of the pitch motion, .295 Hz. 

These steady state oscillations appear quite small until 
it is realized that the maximum roll perturbation that was 
able to be induced during any run was only about . degrees. 
Lower ratios, between induced and background oscillations 
were the norm for both roll and pitch induced motion. The 
background oscillations in pitch and roll therefore could 
not be ignored when analyzing data, and deleterlously 
affected the measurement of decay time constants. 

As stated before the attempt was made to Introduce 
separately pitch and roll motion into the test craft and 
then record the transient response after the forcing action 
had passed or was stopped. After studying the processed 
data it was evident that it had been nearly impossible to 
affect transients in one variable (either pitch or roll) 
without introducing at least some oscillations In^ the other 
variable. In Introducing pitch oscillations roll motion was 
greatly affected. This was caused because the wave front 
caused by the chase boat's wake always struck the test craft 
at an oblique angle and not perpendicular to the axis of 
symmetry of the craft. Only during those runs where the 
XR-3 was sallied was some modicum of success achieved in 
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Isolating roll transient motion. The one run included in 
this report where sallying was used to induce roll transients 
is the one conducted at 11.79 knots. During the other tv;o 
runs included (15.93 and 19.^5 knots) the v;ake crossing 
method v;as used to force oscillations. This resulted In 
runs in which both roll and pitch transients could be 
studied at the same time. 

In making comparison between the recorded data and that 
of the computer simulations both the original XR-3 model as 
presented in Ref. 3 and the new model which incorporates 
the new fan maps, new bow and stern seals and a water slope 
correction were used. This last modification is one which 
has been incorporated into the XR-3 program by LCDR G. T. 
Forbes, USN during his concurrent study and is reported in 
Ref. 5. The water slope modification is used to correct 
the water level in the plenum chamber based upon the fact 
that increased pressure in the plenum will tend to displace 
water after it comes in contact with the bow seal and passes 
into the chamber. This displacement since it does not take 
place instantaneously is a function of both time and pressure 
and hence results in a "slope" in the water as it passes 
through the plenum chamber. The purpose of using both the' 
original and the new models was to gain an appreciation of the 
effects of the modeling changes made. 

During all the runs made at Lake San Antonio with the 
XR-3 the v/ater conditions were fairly calm, but the steady 
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state roll and pitch oscillations were always present In 
the processed data as reported previously. Figures 19 
through 35 which are the graphs for all the runs clearly 
show this oscillation In both pitch and roll. Again It Is 
emphasized that these oscillations precluded exact measure- 
ments of the decay time constants. 

The method used to Induce roll and pitch transient 
oscillations In the computer models was to offset the Initial 
roll and pitch angle from the known steady state angles 
when the program was executed. These offsets would then 
decay as the craft returned to the steady state conditions 
for that speed. The magnitude of these offsets were taken 
from the recorded data. The point in time where this data 
was taken is marked on each graph of the test craft's 
recorded data by a dashed vertical line. 

Prom known loading conditions the craft's weight during 
these runs was 605 O lbs. Its center of gravity was located 
10.08 feet forward of the transOm. Average plenum pressure 
was not available as this data was not recorded. The plenum 
pressure was assumed to be the nominal 2^1.86 PSF as reported 
In Re f . 3 ♦ 

These inputs in the case of the original [Ref. 3] XR-3 
model caused the craft to have a negative steady state trim 
angle which Is readily apparent when looking at the graphs 
of the simulation runs of the original model (see Figures 
22 through 28 for simulation runs at 15.93 knots). 
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Although It was recognized that this negative trim angle was 
not realistic the initial conditions were not varied to get 
a positive trim angle because it was desired to test each 
model with exactly the same physical characteristics. The 
actual value of the steady state trim was not the quantity 
that was to be studied. The quantity that was of interest 
was the transient response of the model as it arrived at 
that steady state value after an initial offset. 

B. RUN DATA ANALYSIS 

The first run v;as a calm water transient roll motion 
study at 11.79 knots. In this run the pitch offset was 
considered negligible. The roll induced oscillation was 
measured to have a negative magnitude of approximately 
.5 degrees. This roll offset was entered into both models 
with the following results. The decay time constant for 
roll motion was ^.6 seconds for the original XR-3 model and 
5.6 seconds for the new model. The oscillation frequency 
was .752 Hz for the original XR-3 and .71^ Hz for the new. 
Measured data gave the actual boat a frequency of roll 
oscillation of .29 Hz and a decay time constant of less 
than 1 second. 

The second run was at 15.93 knots in calm water. This 
run had both roll and pitch perturbations. They were 
minus .25 degrees for pitch and a negative .27 degrees for 
roll. V/lth these offsets the new XR-3 model gave a roll 
time constant of ^.1 seconds and a. pitch time constant of 






3.4 seconds. The frequencies of oscillation for the new 
model were .72 Hz for roll and .44 Hz for pitch. These 
same quantities were measured on the old model to be 

3.1 seconds and 5.0 seconds for the roll and pitch time 
constants respectively, and in order, .76 Hz and .50 Hz for 
the roll and pitch frequencies. The same measured quantities 
from the recorded data gave something less than one second 
for both the roll and pitch time constants and .25 Hz and 
.22 Hz for the oscillation frequencies for roll and pitch. 

The last calm water data was taken at 19.45 knots. The 
measured pitch offset was a positive .12 degree and a 
positive .45 degree for the roll offset. Again these 
values were Imposed upon the initial steady state trim 
conditions to achieve a transient response. 

Prom the two computer simulations the following results 
were obtained. Original model roll frequency — .76 Hz, 
roll decay time constant — 2.7 seconds, pitch frequency — 

.50 Hz, pitch time constant — 4.6 seconds. New model roll 
frequency is .72 Hz and decay time constant 4.3 seconds. 

Pitch frequency is .42 Hz and decay time constant is 

4.4 seconds. The same data from the recorded tapes gave a 
pitch frequency of .28 Hz and with decay time constant of 
less than one second. Roll frequency was measured to be 
.28 Hz, and the roll decay time constant was estimated at 
less than one second. 
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Prom the above data taken on the three runs and the 
computer simulations thereof It Is evident that the computer 
model neither In the original form nor in its new form 
accurately predicts roll or pitch transient motion under 
the conditions studied. Decay time constants for pitch 
and roll were generally four' times too long in both simula- 
tion models. Pitch and roll frequencies of the computer 
models were generally too high by a factor of two. Table 
III gives this transient behavior of the XR-3 and both 
simulation models in tabular form. 

Figure 19 is the actual roll response data for Run 1 
collected at 11.79 knots. Figures 20 and 21 are the 
computer simulation runs for the original and new model 
for this run depicting roll response. Recorded pitch and 
roll motion for Run 2 (15.93 knots) is given in Figures 22, 
23 and 2 ^ and the original model computer simulation of 
roll and pitch is in Figures 25 and 26 respectively. The 
response of the new model for Run 2 is recorded in graphs 
27 and 28. Recorded data for Run 3 is given in Figures 29 
through 31 with the computer model's simulations of pitch 
and roll given in Figures 32 through 35. 

C. SIMULATION STUDIES OF THE TWO MODELS 

In an effort to gain better knowledge of the total 
effect of the new modeling, both the old and the new models 
were given a positive three degree pitch offset from steady 
state calm water conditions. Both models were set at the 



TABLE III. CALM WATER FORCED OSCILLATIONS RESPONSE 



Run 1, Speed 11.79 knots 



Variable 



Measured 

XR3 



Simulated 

Original Model New Model 



Roll Frequency .29 Hz 

Roll decay time constant <1 sec 
Pitch Frequency 
Pitch decay time constant 



.752 Hz .71^ Hz 

4.6 sec 5.6 sec 



Run 2, Speed 15.93 knots 

rfeasured Simulated 



Variable 


XR3 


Original IVfodel 


New Model 


Roll Frequency 


.25 Hz 


.76 Hz 


.72 Hz 


Roll decay time constant 


<1 sec 


3.1 sec 


4.1 sec 


Pitdi Frequency 


.22 Hz 


.50 Hz 


.44 Hz 


Pitch decay time constant 


<1 sec 


5.0 sec 


3.4 



Run 3, Speed 19-45 knots 

ffeasured Simulated 



Variable 


XR3 


Original Model 


New Model 


Roll Frequency 


.28 Hz 


.76 Hz 


.72 Hz 


Roll decay time constant 


<1 sec 


2.7 sec 


4.3 sec 


Pitch Frequency 


.28 Hz 


.50 Hz 


.42 Hz 


Pitch decay time constant 


<1 sec 


4.6 sec 


4.4 sec 
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FIGURE 19 

Roll Transient verses Time> Velocity 11.79 knots 
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FIGURE 20 

Computer Simulation Roll Angle verses Time, 11,79 knots 
Original XR-3 
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FIGURE 21 

Computer Simulation Roll Angle verses Time^ 11.79 knots 
New XR-3 
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FIGURE 22 

Pitch verses Time^ Velocity 15.93 knots 
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FIGURE 23' 

Pitch Transient verses Time^ Velocity 15.93 knots 



82 



w 









FIGURE 2 H 

Roll Transient verses Time^ Velocity 15.93 knots 
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FIGURE 25 

Computer Simulation Pitch Angle verses Time^ 15,93 knots 
Original XR-3 
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FIGURE 26 ’ 

Computer Simulation Roll Angle verses Time^ 15.93 knots 
Original XR-3 
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FIGURE 27 

Computer Simulation Pitch Angle verses Time> 15.93 knots 
New XR-3 
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FIGURE 28 

Computer Simulation Roll Angle verses Time^15.93 knots 
New XR-3 
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FIGURE 29 

Roll Transient verses Time^ Velocity 19.45 knots 
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FIGURE 30 

Pitch verses Time^ Velocity 19.45 knots 
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FIGURE 31 



Pitch Transient verses Time. Velocity 19.^15 knots 
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FIGURE 32 

Computer Simulation Roll Angle verses Time,' 19. A5 knots 
Original XR-3 
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FIGURE 33 

Computer Simulation Pitch Angle verses Time^ 19.45 knots 
Original XR-3 
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FIGURE 3^ 

Computer Simulation Roll Angle verses Time^ 19. A5 knots 
New XR-3 
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FIGURE 35 

Computer Simulation Pitch Angle verses Time^ 19. A5 knots 
New XR-3 
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same Initial conditions that they were during the 19.^5 knot 
run except for the offset. In the case of the new model 
this raised the bow approximately six Inches and since the 
steady state condition gave the bow a depth of about three 
inches with the positive steady state pitch angle this caused 
air gaps to appear beneath the bow seal. In the case of the 
old model where the steady state pitch angle was negative 
no ,alr gaps were created. Figures 36 through 38 present the 
original model’s response and Figures 39 through 4l present 
the new model’s response. 

As can be seen in Figure i|0, the air gap on the new model 
caused wide fluctuations in the plenum pressure during the 
first half-second. This loss of pressure permitted the boat 
to rapidly sink as shown in Figure iJl until the plenum 
pressure built up. When the air gaps closed and the boat 
continued to settle the plenum pressure rapidly rose as the 
plenum volume decreased. The boat exponentially returned 
to steady state draft levels as the pitch motion shown in 
Figure 39 died down. The pitch motion of the new model 
shows the effect of the new seals. The downward movement 
of the bow as seen in Figure 39 is severely limited as the 
bov7 seal wetted length becomes greater and the effect of 
the rigid body analysis comes into play. The first negative 
overshoot is limited to 3^ percent of the original effect. 

The initial decay time constant is 2.5 seconds and the 
pitch frequency is .63 Hz v.'hlch because of the nonlinearity 
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Computer Simulation 3 Degree Pitch Step Response^ 19.^5 knots 
Original XR-3 Pitch verses Time 
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FIGURE 37 

Computer Simulation 3 Degree Pitch Step Response, 19.45 knots 
Original XR-3 Plenum Pressure verses. Time 
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FIGURE 38 

Computer Simulation 3 Degree Pitch Step Response^ 19.^5 knots 
Original XR-3 Draft at C.G. verses Time 
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New XR-3 Pitch verses Time 
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FIGURE HO 

Computer Simulation 3 Degree Pitch Step Response^ 19.^5 knots 
New XR'-3 Plenum Pressure verses Time. 
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FIGURE m 

Computer Simulation 3 Degree Pitch Step Response^ 19.^5 knots 
New XR-3 Draft at C.G. verses Time 
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of the bow and stern seal models changes to a frequency of 
. i|2 Hz at the end of ten seconds. The time constant has 
lengthened to 4.3 seconds in the same period. 

The behavior of the old model is different. Because 
of the high capacity of the fan maps the damping effect of 
plenum pressure is minimal. The draft movement shown in 
Figure 38 is therefore highly underdamped whereas the draft 
motion for the new model did not display any oscillatory 
tendencies (see Figure 4l) and displayed exponential decay 
type response when returning to steady state conditions. 

In pitch motion shown in Figure 36 the original model 
displays a first overshoot magnitude of 78 percent of 
original offset. The frequency of oscillation is . 5 Hz 
and the decay time constant is 4.78 seconds. Both the time 
constant and pitch oscillation frequency remain constant 
throughout the 10 second simulation. 

To further evaluate the differences between the original 
and new models a second simulation run was made at 19.45 
knots. This run was in ahead sea state 2 irregular seas. 
Figures 42 and 43 show that the plenum pressure fluctuations 
in the new model were 6 OO percent greater than in the 
original model. Whereas the plenum pressure varied from 
5 to about 45 PSP in the new model the pressure only varied 
from about 21 to 27 PSP in the original. Draft fluctuations 
were about 300 percent greater in the new model than in the 
original (see Figures 44 and 45). Pitch angle fluctuations 
vrere comparable in both models (Figures 46 and 47). 
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FIGURE l\2 

Computer Simulation Plenum Pressure verses Time> 19. ^^5 knots 
New XR-3^ Sea State 2 
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FIGURE H3 

Computer Simulation Plenum Pressure verses Time^ 19.^5 knots 
Original XR-3^ Sea State 2 
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FIGURE. 

Computer Simulation Draft at C.6. verses Time^ 19,^5 knots 
New XR-3, Sea State 2 
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FIGURE i»5 

Computer Simulation Draft at C.G, verses Time> 19,^5 knots 
Original XR-3, Sea State 2 
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FIGURE il6 • 

Computer Simulation Pitch Angle verses Time^ 19.^5 knots 
New XR-3^ Sea State 2 
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Computer Simulation Pitch Angle verses Time, 19.^5 knots ' 
Original XR-3, Sea State 2 



108 



V. CONCLUSIONS AND RECOMMENDATIONS 



Making accurate measurements from the collected data 
was extremely difficult because of the high relative level 
of steady state oscillation present in both the pitch and 
roll motion. Signal noise coupled with the fact that the 
decay time constants and the periods of the oscillations 
were of the same general magnitude made accurate measurement 
of the decay time constants practically impossible. 

Although accurate measurement of the time constants was 
not possible general observations could be made. It was 
clearly evident that the pitch and roll transients Induced 
in the XR-3 test craft decayed very rapidly and not with 
the highly underdamped characteristics shown by the computer 
simulation models. The measured tape data confirmed the 
observations made aboard the XR-3 while conducting these 
runs that the craft come to steady state conditions with 
only one or two overshoots after the forcing action had 
ceased. Roll motion transients were normally completed in 
less than tv;o cycles. Both the longitudinal and athwartshlps 
stability were highly speed dependent gaining stability as 
speed Increased. 

These observations were not based solely on "feel" but 
were made by sighting across the deck of the XR-3 and 
observing the motion of the deck edge relative to the shore 
line . 
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simulation runs have shown that for large perturbations 
the transient response of the new model Is quite different 
from the original. Since it was not possible to excite such 
large perturbations in the actual test craft it remains to 
be verified that the new bow and stern seal models more 
accurately predict the forces generated by the actual seals 
than does the original seal model. 

The new seal fan maps are more accurate than the original 
since they are based on actual verified manufacturer's data. 
Leakage areas are now proportional to those of the actual 
craft . 

The fact that the present computer models of the XR-3 
both old and new do not appear to accurately represent the 
actual transient behavior of the XR-3 in pitch and roll 
opens many areas which should be investigated for possible 
errors. As both the pitch and roll transients appear to be 
in error the sidewall subroutine becomes a prime area of 
investigation as it contributes major components to both 
the roll and pitch force moments. 

Figure *l8 shows one area which should be explored and 
that is the area of hydrodynamic lift of the sidewalls. 

Figure is a simplified drawing of the underside of 
the XR-3. The flat sections on the bottom of each 
sidev;all are cross-hatched for clarity. At high 
velocity (20-25 knots) these areas become effective 
planing surfaces. In this program the lift of the 
sidewall sections due to planing does not appear to take 
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into account the underwater body shape. This shape changes 
radically from bow to stern. Whether or not this shape has 
significant effect on roll and pitch transient behavior is 
not known but should be investigated. Another aspect of 
the sidev/all modeling which should be investigated is the 
inclusion of forces due to the displacement of water per 
unit of time. A term such as this would not only have terms 
based on draft and pitch angle but would also include terms 
which took into account the rate of change of displacement 
and the rate of change of the pitch angle. The addition of 
first derivative terms such as these would have a definite 
effect on the transient behavior of the sidewalls. 

A third area of study concerning the sidewalls is the 
effect of sidewall deadrise on turning performance. Modeling 
of the deadrise forces would be incorporated in the cross- 
flow force modeling. This deadrise could have significant 
effect on the outboard sidewall during a turn determining 
whether it has a tendency to rise and ride over the water 
as the craft skids or dig in and cause the craft to undergo 
violent roll motion as reported in Ref. 6. 

Another area of study which needs investigation is the 
effect of plenum pressure on the transient response (of the 
computer model). If it can be shown that the computer 
model as presently configured can be made to duplicate 
actual measured craft transient motion by just adjusting 
the plenum pressure, further modifications of the program 
m.ay be unnecessary. Further transient motion studies are 
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needed where both the plenum pressure and the seal pressures 
are included in the recorded data. Because of the stochastic 
nature of the actual measured data many duplicate runs will 
be needed to arrive at accurate averages of the test craft's 
oscillation frequencies and the decay time constants. These 
averages should then be used in comparing the actual craft 
with the computer model. 

As can be seen in the recorded data presented herein a 
fruitful area of study would be in reducing the signal noise 
on the recorded data. This could be in the recording phase 
or processing phase of the data. Thorough tests under 
controlled conditions need to be made of the actual transfer 
characteristics of the different sensors installed in the 
test craft. 

Earlier in this report it was mentioned that there 
existed steady state roll and pitch oscillations while the 
boat was traveling straight on calm water. These oscillations 
could exist and could be the result of the craft being in 
some limit cycle type of motion because of the inherent 
nonlinearity of vehicles operating in fluid mediums. On 
the other hand, it could be the result of instrumentation 
for example if the vertical seeking gyro from which all roll 
and pitch motions are derived was precessing about the true 
vertical then it could impose a sinusoidal oscillation 
signal on both the pitch and roll sensors that would show 
up in the recorded data. 
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•lilllllllinsm 









In conclusion it must be stated that there are still 
many areas in this subject of sidewall modeling of the 
XR-3 which should and need to be explored before the XR-3 
model as originally stated in Ref. 3 and modified herein 
can be accepted as a reasonably accurate representation of 
the actual XR-3 test craft. Further investigation into 
the data collection and processing is also needed to refine 
the processed data to the point where it is accurate 
enough for dynamic studies. 
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APPENDIX A 



LISTING OF BOWSEAL SUBROUTINE 'bOWSl' 



The entering arguments needed for the new bow seal 
routine are the same as presented In Ref. 3 for the original 
bow seal except that Instead of entering the angle between 
the leading edge and craft vertical In columns 46-55 j the 
length of the rear cables Is entered (In feet). An additional 
card Is needed to enter the length of the middle support 



cables . 


The new 


data card should 


be entered as shown below. 


CARD 


COLUMNS 


FORMAT 


ENTRY 


1 


1- 5 


I 


00601 — control tag 




6-15 


F 


X coordinate of seal 
hinge (ft) forward 
of transom 




16-25 


F 


Seal leakage orifice 
coefficient 




26-35 


F 


Pressure differential 
between bow seal and 
plenum chamber (PSF) 




36-45 


F 


Z coordinate of seal 
hinge (ft) above keel 




46-55 




Length of rear support 
cable (ft) 

( 1.875 Is full down) 




56-65 




Length of the seal 
leading edge (ft) 




66-75 




2 

Base leakage area (ft ) 
Nominal .08 ft^) 


2 


1- 5 




00602 — control tag 




6-15 




Length of middle support 



cable (ft) 

(1.1 B?5 Is full down) 
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APPENDIX B 



LISTING OF STERN SEAL SUBROUTINE STNSL 



The entries for this routine are the same as the original 
stern seal except that in columns ^6-55 the rear support 
cable length is entered Instead of the seal angle. 

Block 5 entries are as follows: 



CARD 


COLUMNS 


FORMAT 


ENTRY 


1 


1- 5 


I 


00500 — control tag 




6-15 


F 


X coordinate of seal hinge 
(ft) forward of transom 




16-25 


F 


Z coordinate of seal hinge 
above keel (ft) 




26-35 


F 


Base leakage area (ft ) 
3.79 nominal 




36-il5 


F 


Seal leakage orifice 
coefficient 




il6-55 


F 


Length of rear support 
cable (ft) 1.875 is full 
down 




56-65 


F 


Pressure differential 
between stern seal and 
plenum chamber (PSF) 




66-75 


F 


Length of the leading edge 
of the seal (ft) 
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APPENDIX C 



LISTING OP PAN SUBROUTINE PAN' 

The entries for the fan subroutine remain the same. 
The nominal fan RPM Is 8000 RPM vice 2615 RPM. Printed 
below are copies of the appropriate entries for the bow 
seal 3 stern seal and main plenum fans. 



01901 2.0 
- 10.406 
23.414 
52.030 
22.7 
12.20 
-.020 

01902 2.0 
- 10.406 
23.414 
52.030 
22.7 
12.20 
-.020 

01903 1.0 
- 10.406 
23.414 
52.030 
25 . 70 
16.40 
0.00 



8000 


.0 19.0 


1 


- 5.203 


0.0 


5.203 


26.015 


28.617 


31.218 


57.233 


62.436 




22.067 


20.67 


19.60 


10.00 


7.80 


6 . 40 


-1 .06 


- 2.03 




8000 


1.0 19.0 


1 


- 5.203 


0.0 


5.203 


26.015 


28.617 


31.218 


57.233 


62.436 




22.067 


20.67 


19.60 


10.00 


7.80 


6 . 40 


- 1.06 


- 2.03 





- 5.203 

26.015 

57.233 

24.33 

12.60 

- 1.06 



BOOO.O 



19.0 



0.0 
28.617 
62.436 
23.167 
10.367 
- 2.03 



1.0 

5 . 203 
31.218 

21.83 

8.033 



10.406 

33.820 

18.367 

5.20 



10.406 

33.820 

18.367 

5.20 



10.406 

33.320 

20.61 

6.331 



15.609 

36.421 

17.133 

4.33 



15.609 

36.421 

17.133 

4.33 



15.609 

36.421 

19.33 

5.073 



18.211 

41.624 

16.467 

2.60 



18.211 

41.624 

16.467 

2.60 



18.211 

41.624 

13.85 

3.035 



20.812 

46.827 

15 . 10 
.983 



20. 812 
46. 827 

15.10 

.983 



20.812 

46.827 

17.73 

1.60 
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APPENDIX D 



LISTING OF INITIAL CONDITION SUBROUTINE INCOn' 



The INCON subroutine was modified to permit it to 
calculate the wave components of the various sea states if 
desired instead of having to enter each component Individually. 
Block 11 options 1 and 2 remain identical. New options 



3 and 
BLOCK 


4 have been 
11 - OPTION 


added and 
TAG 3 


CARD 


COLUMNS 


FORMAT 


1 


1- 5 


1 




6-15 


F 




16-25 


F 




26-35 


F 




36-45 


F 




46-55 


F 



BLOCK 


11 - OPTION 


TAG 4 


CARD 


COLUMNS 


FORMAT 


1 


1- 5 


1 




6-15 


F 




16-25 


F 




26-35 


F 



are entered as follows: 



ENTRY 

01103 control and option tag 

Number of wave components 
maximum of 10 

Initial wave heading (deg) 
(l80 deg = ahead waves) 

Average height of 1/3 
highest waves (ft) 

Shortest significant period 
(sec) 

Longest significant period 
(sec) 



ENTRY 

01104 control and option tag 

Niimber of wave components 
maximum of 10 

Initial wave heading (deg) 
(l 80 deg = ahead v;aves) 

Average height of 1/3 
highest vjaves (ft) 
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CARD 


COLUMNS 


FORMAT 


ENTRY 






1 


36-H5 


P 


Lowest 

(Hz) 


significant 


frequency 




H 6-55 


P 


Highest 

(Hz) 


significant 


frequency 
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LISTING OF INITIAL CONDITION SUBROUTINE INCON 
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